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I. INTRODUCTION 



The determination of the properties of hadronic matter at finite temperature and den- 
sity is a fundamental problem in nuclear physics. In particular, the study of the liquid-gas 
phase transition in medium energy heavy-ion collisions is of considerable interest. Many 
intermediate-energy collision experiments have been performed [1] to investigate the un- 
known features of the highly excited (hot) nuclei formed in such collisions [2,3]. Theoret- 
ically, much effort has been devoted to studying the equation of state for nuclear matter 
and to discussing the critical temperature, T^. Recently, Natowitz et al. obtained the lim- 
iting temperature by using a number of different experimental measurements [4]. From 
these observations the authors extracted the critical temperature of infinite nuclear matter, 
Tc = 16.6 ± 0.86 MeV [5]. We can expect that further experiments may eventually yield the 
limiting temperature of hypernuclei and the critical temperature for infinite strange hadronic 
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matter. It is therefore interesting to study the hquid-gas phase transition of strange hadronic 

matter theoretically. 

Exploring systems with strangeness, especially with large strangeness fraction, has at- 
tracted a lot of interest in recent years. Such a system has many astrophysical and cosmolog- 
ical implications and is indeed interesting by itself. There are many theoretical discussions 
for both strange hadronic matter [6,7] and strange quark matter [8]- [11]. However, most 
discussions are at zero temperature. The properties of strange hadronic matter at finite 
temperature have not been studied very much yet. Unlike symmetric nuclear matter, for 
strange hadronic matter, there are two conserved charges, baryon number and strangeness. 
Glendenning [12] first discussed the phase transition with more than one conserved charge 
in general and applied it to the possible transition to quark matter in the core of neutron 
stars. Miiller and Serot [13] discussed asymmetric nuclear matter, which has two conserved 
charges (baryon number and isospin), using the stability conditions on the free energy, the 
conservation laws and the Gibbs criterion for the liquid-gas phase transition. The liquid-gas 
phase transition of asymmetric nuclear matter was also discussed in effective chiral mod- 
els in Refs. [14,15]. It was found that the critical temperature decreases with increasing 
asymmetry parameter a. 

For strange hadronic matter, the method is similar to that for asymmetric nuclear matter. 
In both cases there are two conserved charges. Recently, Yang et al. [16] used the extended 
Furnstahl-Serot-Tang (FST) model [17] to discuss the hquid-gas phase transition of strange 
hadronic matter. The original FST model satisfies the SU (2) chiral symmetry and was first 
applied to study nuclear matter [18]. In the extended FST model [16], the author found 
a so called critical pressure above which the liquid-gas phase transition cannot exist. As a 
result no critical strangeness fraction was obtained for a given temperature. This critical 
pressure exists in finite nuclei together with a limiting temperature, T^^. However, in infinite 
hadronic matter, physically, there should be no such critical pressure, since the pressure of 
the system can be much higher than the critical pressure. For strange matter with fs = 2, 
i.e. pure S matter, from Fig. 1 of Ref. [16] one sees that the pressure does not increase 
monotonically with density. This means that at this temperature, pure S matter can be in 
hquid-gas phase coexistence. Therefore, the binodal p — /jl diagram at temperature T = 10 
MeV should terminate at = 2. We will reconsider this problem and show how the critical 
strangeness fraction can be obtained. 

To study the properties of hadronic matter, we need phenomenological models since 
QCD cannot yet be used directly. The symmetries of QCD can be used to constrain the 
hadronic interactions and models based on SU (2) ^ x SU (2)^^ symmetry and scale invariance 
have been proposed. These effective models have been widely used to investigate nuclear 
matter and finite nuclei, both at zero and at finite temperature [18]- [20]. Papazoglou et al. 
extended the chiral effective models to SU (3)l x SU (3)r, including the baryon octets [21,22]. 
As well as models based on hadronic degrees of freedom, there are some other models based 
on quark degrees of freedom, such as the quark meson couphng model [23,24], the cloudy 
bag model [25], the NJL model [26] and the quark mean field model [27], etc.. Recently, we 
proposed a chiral SU{3) quark mean field model and investigated the properties of hadronic 
matter as well as quark matter [28]- [32]. This model is quite successful in describing the 
properties of nuclear matter [28], strange matter [30,31], finite nuclei and hypernuclei [32] 
at zero temperature. In this paper, we will apply the chiral SU (3) quark mean field model 
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to finite temperature and study the liquid-gas phase transition of strange hadronic matter. 

The paper is organized as follows. The model is introduced in section II. In section III we 
use the model to investigate strange hadronic matter at finite temperature. The numerical 
results are discussed in section IV and section V summarises our findings. 



II. THE MODEL 

Our considerations are based on the chiral SU(3) quark mean field model (for details 
see Refs. [30,32]), which contains quarks and mesons as basic degrees of freedom. Quarks 
are confined in baryons by an effective potential. The quark meson interaction and meson 
self-interaction are based on SU (3) chiral symmetry. Through the mechanism of sponta- 
neous chiral symmetry breaking, the resulting constituent quarks and mesons (except for 
the pseudoscalars) obtain masses. The introduction of an explicit symmetry breaking term 
in the meson self-interaction generates the masses of the pseudoscalar mesons which satisfy 
partially conserved axial-vector current (PCAC) relations. The explicit symmetry break- 
ing term of the quark meson interaction leads in turn to reasonable hyperon potentials in 
hadronic matter. For completeness, we introduce the main concepts of the model in this 
section. 

In the chiral limit, the quark field q can be split into left and right-handed parts Ql and 
Qr- Q = Ql + Qr- Under SU{3)l x SU{3)r they transform as 

q'l = Lql, q'n = RqR. (1) 
The spin-0 mesons are written in the compact form 

1 ^ 

M(M+) = E±m= ^^(a"±i7r'^)A", (2) 

v2 „=o 

where a" and tt" are the nonets of scalar and pseudoscalar mesons, respectively. A" (a = 
1, ...,8) are the Gell-Mann matrices, and A° — \/|-^- The alternative plus and minus signs 
correspond to M and M+. Under chiral SU{3) transformations, M and transform as 
M ^ M' = LMR+ and M+ M+' = RM+L+. In a similar way, the spin-1 mesons are 
introduced through: 

Ur,) = ^ (^M ± \) = ^ E (< ± <) A" (3) 

with the transformation properties: l'^ = Ll^L~^, r'^ = Rr^R^. The matrices S, 

n, Vn and can be written in a form where the physical states are explicit. For the scalar 
and vector nonets, we have the expressions 



1 s 
E = ^ ^ a" A" 
v2 „=o 



K*- K*° C I 



(4) 
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a=0 



pt ^ 



(5) 



Pseudoscalar and pseudovector nonet mesons can be written in a similar fashion. 
The total effective Lagrangian has the form: 

^eff = ^qO + ^qM + ^SS + ^VV + ^xSB + -^Am^ + ^hi + ^ci 



(6) 



where CqQ — q i'j'^d^ q is the free part for massless quarks. The quark-meson interaction CgM 
can be written in a chiral SU (3) invariant way as 



\a=0 



a=0 



(7) 



In the mean field approximation, the chiral-invariant scalar meson >Cee and vector meson 
Cvv self-interaction terms are written as [30,32] 



'EE 



1 



X 



-hx - 7% In— + -X In- ^ , 



ixf 

2X^ 



-Cyy = ^ ^ + my + m|0') + ^4 {uJ^ + 60;^' + + 20^) , 



(8) 
(9) 



where S = 6/33; ctq, Co and Xo are the vacuum expectation values of the corresponding mean 
fields (7, ( and x- The Lagrangian jC^sb generates the nonvanishing masses of pseudoscalar 
mesons 



-XSB 



r 

Xo 



(10) 



leading to a nonvanishing divergence of the axial currents which in turn satisfy the relevant 
PCAC relations for tt and K mesons. Pseudoscalar and scalar mesons as well as the dilaton 
field X obtain mass terms by spontaneous breaking of chiral symmetry in the Lagrangian 
(8). The masses of the m, d and s quarks are generated by the vacuum expectation values 
of the two scalar mesons, a and C,. To obtain the correct constituent mass of the strange 
quark, an additional mass term has to be added: 



^Am, = -Am^qSq, 



(11) 



where S — ^ (^I — XsV^^ — diag(0, 0, 1), is the strangeness quark matrix. Through these 
mechanisms, the quark constituent masses are finally given by 



mu = md = -^^0 and rUs = -QsCo 



Am. 



(12) 
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where Qs and Am^ arc chosen to yield the constituent quark mass in vacuum - in our case, 
iriu = tTid = 313 MeV and m.^ = 490 MeV. In order to obtain reasonable hyperon potentials 
in hadronic matter, it has been found necessary to include an additional coupling between 
strange quarks and the scalar mesons a and ( [30] . This term is expressed as 

Ch ^ {hia + h20ss. (13) 

In the quark mean field model, quarks are confined in baryons by the Lagrangian Cc — 
—^Xc^ (with Xc given in Eq. (14), below). The Dirac equation for the quark field ^ij, 
under the additional infiuence of the meson mean fields, is given by 

— * 

-id -V + Xc{r) + (3711* 

where a = 'y^'y , /5 = 7° , the subscripts i and j denote the quark i {i = u, d, s) in a baryon of 
type j (j = A^, A, S, S) and Xc(^) is a confining potential - i.e. a static potential providing 
confinement of quarks by meson mean-field configurations. The quark effective mass, m*, 
and energy e* are defined as 

m* = -gia - giC + rriio (15) 

and 

e* = - 9> - g]t>(t> , (16) 

where is the energy of the quark under the infiuence of the meson mean fields. Here 
rriio = for i = u,d (nonstrange quark) and mjo — ArUs — 29 MeV for i = s (strange 
quark). Using the solution of the Dirac equation (14) for the quark energy e* it has been 
common to define the effective mass of the baryon j through the ansatz: 

M; = ^Ef- < > , (17) 

where E* = Ylii^ij^ + Ejgpin is the baryon energy and < p*cm > i^ ^he subtraction of the 
contribution to the total energy associated with spurious center of mass motion. In the 
expression for the baryon energy riij is the number of quarks with fiavor "i" in a baryon 
with fiavor j, with j = N {p, n} , S {S^*", E°} , S {E^, S~} , A and Ejapin is the correction to 
the baryon energy which is determined from a fit to the data for baryon masses. There is 
an alternative way to remove the spurious c. m. motion and determine the effective baryon 
masses. In Ref . [33] , the removal of the spurious c. m. motion for three quarks moving in a 
confining, relativistic oscillator potential was studied in some detail. It was found that when 
an external scalar potential was applied, the effective mass obtained from the interaction 
Lagrangian could be written as 

M* = Y.n^,el-E% (18) 

i 

where E^ was found to be only very weakly dependent on the external field strength. We 
therefore use Eq. (18), with E^ a constant, independent of the density, which is adjusted to 
give a best fit to the free baryon masses. 



(14) 
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Using the square root ansatz for the effective baryon mass, Eq. (17), the confining po- 
tential Xc is chosen as a combination of scalar (S) and scalar-vector (SV) potentials as in 
Ref. [30]: 

Xc{r) = \[xl{r) + xT{r)] (19) 

with 

xl{r) = -^Kr\ (20) 

and 

xf{r)^\Kr\l + ^'). (21) 

On the other hand, using the linear definition of effective baryon mass, Eq. (18), the confining 
potential Xc is chosen to be the purely scalar potential Xc(^)- The coupling kc is taken as 
kc — 1 (GeV fm~^), which yields baryon mean square charge radii (in the absence of a pion 
cloud [34]) around 0.6 fm. 

The properties of infinite nuclear matter and finite nuclei were calculated with these two 
treatments of effective baryon mass in Ref. [35]. As we have explained there, the linear 
definition of effective baryon mass has been derived using a systematic relativistic approach 
[33] , while to the best of our knowledge no equivalent derivation exists for the square root 
case. For high baryon density, the predictions of these two treatments are quite different. 
Many physical quantities change discontinuously at some critical density in the case of square 
root ansatz, while the linear definition of baryon mass yields continuous behavior for high 
density nuclear matter. Both treatments of the spurious c. m. motion fit the saturation 
properties of nuclear matter and therefore, for densities lower than the saturation density, 
these two treatments give reasonably similar results. In this paper, we will discuss the 
liquid-gas phase transition of strange hadronic matter with both treatments. We prefer the 
linear form because it has been derived. The square root case is reported here because it is 
widely used and in fact produces similar results in some regions. However, where they differ 
we believe that the linear form is the more reliable. 



III. STRANGE HADRONIC MATTER AT FINITE TEMPERATURE 

Based on the previously defined interaction, the Lagrangian density for strange hadronic 
matter is written as 

-g^i'Bit.^BUJ^ - g^^Bi^.i'B4>^ + Cm, (22) 

where 

Ffj,u = dfj,uj„ - d^LJf, and Sf,i, = d^,(p^ - dyCp^. (23) 
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The term Cm represents the interaction between mesons which includes the scalar meson 
self- interaction £sS) the vector meson self-interaction Cyv and the explicit chiral symmetry 
breaking term C^^sbi all defined previously. The Lagrangian includes the scalar mesons cr, C, 
and X, and the vector mesons uj and 0. The interactions between quarks and scalar mesons 
result in the effective baryon masses M%, where subscript B labels the baryon B — N, K,Ti 
or S. The interactions between quarks and vector mesons generate the baryon- vector meson 
interaction. The corresponding vector coupling constants and are baryon dependent 
and satisfy the relevant SU{2i) relationships. In fact, we find the following relations for the 
vector coupling constants: 

9t = 9l = 2gl = lg!: and ^ = = = (24) 

At finite temperature and density, the thermodynamic potential for strange hadronic 
matter is defined as 

^ = - E 1^ r '^'^W + e-^^^W-'^^)/'^-^) (25) 

A— AT \ (-^TT) "^0 I 



j=JV,A,E 



+ ln(l + e-(^^(^)+'^^)/^-^)}-£M, 



where E*{k) = \j M*^ -|- 1^ and Qj is the degeneracy of baryon j {gN,^ = 2, = 1 and g-£ — 

3). The quantity Uj is related to the usual chemical potential iij by uj = fj^j — g^uj — (^^0. The 
energy per unit volume and the pressure of the system can be derived as e — Q — + VjPj 
and p = —fl, where pj is the baryon density. 

The mean field equation for meson 0j is obtained by the formula dfl/d(j)i — 0. For 
example, the equations for a, ( are deduced as: 

^ ^ 3(7 Xo 



m,„u 



XoJ da ,=^s,s 9^ 



+ E ^<^.^. >=0, (26) 



where 



a- foo k'^M* 

< i/;,^',- >= ^ dk^j^[nj(k)+nj(k)]. (28) 

In the above equation, nj{k) and nj{k) are the baryon and antibaryon distributions, respec- 
tively, expressed as 
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exp 



{E;{k) - uj) /ksT 



+ 1 



and 



(29) 



nj{k) 



exp 



{E*{k) + v,) /ksT 



+ 1 



The equations for vector mesons cu and are expressed as 

AO j=Ar,A,S,S 



(30) 



(31) 



^mj0 + 8^4<; 
Xo 



where pj is the density of baryons of type j, expressed as 



9j 



dkk^ [nj{k) — nj{k)] . 



(32) 



(33) 



Let us now discuss the hquid-gas phase transition. For strange hadronic matter, we 
follow the thermodynamic approach of Refs. [12] and [13]. The system will be stable against 
separation into two phases if the free energy of a single phase is lower than the free energy 
in all two-phase configurations. This requirement can be formulated as [13] 



F{T,p)<{l-X)F{T,p') + XF{T,p"), 



(34) 



where 



p = (1 - A)p' + Ap", 0<A<1, 



(35) 



and F is the Helmholtz free energy per unit volume. The two phases are denoted by a 
prime and a double prime. If the stability condition is violated, a system with two phases 
is energetically favorable. The phase coexistence is governed by the Gibbs conditions: 



l,'AT,p')^l,''iT,p") (i = Ar,A,E,S), 



(36) 



p'{T,p')^p"{T,p"), (37) 

where the temperature is the same in the two phases. The chemical potentials of the baryons 
satisfy the following relationship: 

A^A = A*s = {pn + A*s)/2. (38) 

Therefore, there are only two independent chemical potentials for the four kinds of baryons. 
They are determined by the total baryon density, ps and the strangeness fraction, fg, which 
are defined as pb = {pn + Pa + Pe + Ps) and fs = (pa + Ps + '^Ps)/pb- 
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IV. NUMERICAL RESULTS AND DISCUSSIONS 



The parameters in this model were determined by the meson masses in vacuum and the 
properties of nuclear matter which were hsted in table I of Ref. [35]. We now discuss the 
liquid-gas phase transition of strange hadronic matter. In Fig. 1, we plot the pressure of the 
system versus baryon density for various strangeness fractions, fs, at temperature T = 15 
MeV for the square root ansatz of the effective baryon mass (Eq. (17)). For nonstrange 
hadronic matter, the p — Pb isotherms exhibit the form of two phase coexistence with an 
unphysical region. The nuclear matter can be in a state of liquid-gas coexistence at this 
temperature. With increasing fg, the pressure will increase. At a particular value of fs, 
the pressure will increase monotonically with increasing density. As we will see later, the 
the strangeness fraction is different in the liquid and gas phases. Therefore, the system 
can still be in liquid-gas coexistence, even though the pressure increases monotonically with 
density. The unphysical region appears again in the range of 1.0 < < 1.6. It is obvious 
that the behavior of the pressure of strange hadronic matter is not monotonic with fs- For 
the linear definition of effective baryon mass (Eq. (18)), the results are plotted in Fig. 2. 
At small strangeness fraction, say fs < 0.4, there are unphysical regions. In the range 
0.4 < fs < 1.75, the pressure increases monotonically with increasing density, while for 
fs > 1.75, the unphysical regions appear again. 




FIG. 1. The pressure of strange hadronic matter p versus baryon density ps with different 
strangeness fraction fs at temperature T = 15 MeV in the case of square root ansatz of effective 
baryon mass. 

As we pointed out earlier, there are two independent chemical potentials for the baryons. 
We now show how the Gibbs conditions can be satisfied. As an example, we plot the 
chemical potentials of nucleon and A versus fs at temperature T = 15 MeV and pressure 
p — 0.23 MeV-fm~^ with the square root ansatz for the effective baryon mass in Fig. 3 (For 
convenience, we use the reduced chemical potential which is defined Mj). The 

solid and dashed lines are for nucleon and A, respectively. The Gibbs equations (36) and 
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(37) for phase equilibriTim demand equal pressure and chemical potentials for two phases 
with different concentrations. The desired solution can be found by means of the geometrical 
construction shown in Fig. 3, which guarantees the same pressure and chemical potentials of 
nucleon and A in the two phases with different fg. Due to the chemical relationship between 
the baryons, the chemical potentials of E and S are also the same in the two phases. 



> 

<u 




pB(fm- ) 

FIG. 2. The pressure of strange hadronic matter p versus baryon density ps with different 
strangeness fraction fg at temperature T = 15 MeV in the case of linear definition of effective 
baryon mass. 
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FIG. 3. Geometrical construction used to obtain the chemical potentials and strangeness fraction 
in the two-phase coexistence at temperature T = 15 MeV and p = 0.23 MeVfm"^ The solid and 
dashed hues are for nucleon and A, respectively. 
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For asymmetric nuclear matter, there is only one kind of solution which satisfies Eqs. 
(36) and (37) at the given pressure and temperature. For hadronic matter, there is another 
solution of the Gibbs equations with higher pressure at the same temperature. We show 
in Fig. 4 the chemical potentials versus at pressure p — 0.27 MeV-fm~^ with the same 
temperature as in Fig. 3. From the geometrical construction, one can see that the difference 
of fs between the two phases is very small. We can discuss in some detail the reason 
why there is only one solution for asymmetric nuclear matter and two solutions for strange 
hadronic matter. This is because for asymmetric nuclear matter with a = {Pn~ Pp) / {Pn+ Pp) , 
the a-dependence-behaviour of the pressure is monotonic - as can be seen clearly in Fig. 3 
of Ref. [14]. However, in the case of strange hadronic matter, the dependence on fg is 
not monotonic. For example, for temperature T = 15 MeV, the point where dp/ dps = 0, 
appears in the two regions of fs - i.e. < < 0.2 and 1.0 < fg < 1.6. This means that the 
system can be in liquid-gas coexistence phase. 



> 




1.36 1.38 1.40 1.42 1.44 1.46 1.48 1.50 
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s 

FIG. 4. Geometrical construction used to obtain the chemical potentials and strangeness fraction 
in the two-phase coexistence at temperature T = 15 MeV and p = 0.27 MeVfm""^ The solid and 
dashed lines are for nucleon and A, respectively. 

The pairs of solutions shown in Fig. 3 with small strangeness fraction form a binodal 
curve which is plotted in Fig. 5. There is a critical point. A, where the pressure is about 
0.244 MeVfm"^ with the corresponding strangeness fraction fg = 0.24. The binodal curve 
is divided into two branches by the critical point. One branch corresponds to the high 
density (liquid) phase, the other corresponds to the low density (gas) phase. Assume the 
system is initially prepared in the low density (gas) phase with fg = 0.2. When the pressure 
increases to some value, the two-phase region is encountered at point a and a liquid phase 
at b with a low fg begins to emerge. As the system is compressed, the gas phase evolves 
from point a to c, while the liquid phase evolves from b to d. If the pressure of the system 
continues to increase, the system will leave the two-phase region at point d. The gas phase 
disappears and the system is entirely in the liquid phase. This kind of phase transition is 
different from the normal first order phase transition where the pressure remains constant 
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during phase transition. If the strangeness is larger than 0.24, there is no phase transition 
between hquid and gas phases. Therefore, for a given temperature there exists a critical 
strangeness fraction, above which the system can only be in the gas phase. In other words, 
for a system with a fixed strangeness fraction fs there exists a critical temperature, above 
which the system cannot change completely into liquid phase however large the pressure. 



> 



0.25 



0.24- 



0.23- 



0.22- 



0.21 -■ 




0.00 



0.30 



FIG. 5. The first binodal curve with smaller strangeness fraction at temperature T 
The points a through d denote the liquid-gas phase transition. A is the critical point. 



15 MeV. 



The solutions shown Fig. 4 with higher strangeness fraction form another binodal curve 
and we plot it in Fig. 6. As in Fig. 5, there are two branches divided by the critical points 
B and C. One branch is for liquid phase , the other is for gas phase . One can see that the 
difference of fg between liquid and gas phase is very small. At the critical points B and C 
the strangeness fraction is about 1.04 and 1.66, respectively. If fg is smaller than 1.04 but 
larger than 0.24, or fg > 1.66, the system can only be in the gas phase. At the critical points 
B and C, the strangeness fraction of hquid and gas phase is the same. The liquid-gas phase 
transitions at these two points are the same as the symmetric nuclear matter. The pressure 
maintains the constant, fs does not change during the phase transition. 
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For the linear case, there are also two kinds of solutions which satisfy the Gibbs equations 
at T = 15 MeV. However, the third critical point C disappears because it moves to the 
point with fg = 2. This is because at this temperature, the system can be in the two phase 
coexistence when fg > 1.75. This can be seen clearly from Fig. 2 where the unphysical 
region exists for very large strangeness fraction fg. 

Therefore, for the case of square root ansatz, altogether there are three critical 
strangeness fractions for strange hadronic matter at temperature T = 15 MeV, while for the 
linear definition of the effective baryon mass, there are only two critical points. In Fig. 7, 
we plot the critical temperature versus strangeness fraction. For the square root case, the 
critical temperature first decreases with increasing fg. When 0.65 < fs < 1-3, increases 
with fs- For fg greater than 1.3, Tc decreases with fg again. In the range, 14.6 < T < 15.3 
MeV, there are three critical strangeness fractions for a given temperature. For example, at 
T = 15 MeV, the three critical values of /, are about 0.24, 1.06 and 1.66. If 0.24 < /, < 1.04, 
the strange hadronic matter can only be in the gas phase. When 1.04 < fg < 1.66, the sys- 
tem can be in a state of liquid-gas coexistence, while for fg > 1.66, the system can once 
again only be in the gas phase. The critical temperature for strange hadronic matter with 
/s = 2 is about 13.6 MeV. For the case of the linear definition of effective baryon mass, 
decreases with increasing fg if fg is smaller than 1.1. If is larger than 1.1, increases 
with fg till fg — 2.0. Compared with the square root case, Tc changes more quickly with fg 
in this case. The highest and lowest Tc is 17.9 and 12.0 MeV with the corresponding fg — 
and fg — 1.1. 
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FIG. 7. The critical temperature Tg versus strangeness fraction fg- The sohd and dashed hues 
correspond to the square root and hnear case, respectively. 



V. SUMMARY 

In this paper, we applied the chiral SU{3) quark mean field model to investigate the 
properties of strange hadronic matter at finite temperature and density. All the parameters 
have been determined in earlier papers and there is no further parameter to be adjusted. 
The model works very well at zero temperature. The saturation properties and compression 
modulus of nuclear matter are reasonable. The hyperon potentials are close to the empirical 
values for strange hadronic matter. The results of finite nuclei and hypernuclei are also 
consistent with the experiments. In this paper, the liquid-gas phase transition of strange 
hadronic matter was studied in this model. For strange hadronic matter, there are two 
independent conserved charges. The system will be in the liquid-gas phase coexistence if 
the pressure and the chemical potentials of all the baryons are the same in the two phases. 
During the phase transition, the strangeness fraction fg of liquid and gas phases is different 
and changes during phase transition, though the total fg is conserved. We found that 
there are two branches of solutions which satisfy the Gibbs equations at some range of 
temperature. One corresponds the phase transition at small fs, while the other corresponds 
the phase transition at large fg. For the square root ansatz of effective baryon mass, there 
are three critical strangeness fractions during 14.6 < T < 15.3 McV. For the linear definition 
of effective mass, there are two critical points when 12.0 < T < 15.9 MeV since the third 
one disappears and moves to the point with fg = 2. The difference of fg in the two phases 
is small, especially in the case of higher strangeness fraction case. The critical temperature 
Tc does not change monotonically with fg. For the square root case, if fg < 0.65 or fg > 1.3, 
Tc decreases, while for 0.65 < fg < 1.3, Tc increases with increasing fg. For the linear case, 
Tc first decreases and then increases with increasing fg. The minimum critical temperature 
is about 12.0 MeV with /, = 1.1. 
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